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We describe the interaction stimulated relaxation in the ensemble of two- 
level systems, responsible for low temperature kinetics and thermodynamics 
properties of amorphous solids. This relaxation gets significant at sufficiently 
low temperature when phonons are substantially frozen out. We show that in 
the realistic experimental situation the measuring field strongly accelerates 
the interaction stimulated relaxation. The characteristic temperature and 
field dependences of the relaxation rate are found when the rate is affected 
both by the interaction between two level systems and by the external field. 
PACS numbers: 61.43.Fs, 75.50.Lk, 11.22. Ch 

1. Introduction 

For a long time the standard model of non-interacting two-level-systems 
(TLS) 1 has served as a good background for understanding experimental 
data in glasses at low temperatures T < IK. The further investigations 
revealed that below lOOmK almost all measurements in dielectric glasses 
dealing with their relaxation properties 2-6 cannot be treated ignoring the 
TLS interaction. For this reason one can suppose that that below WOmK 
TLS manifest the collective behavior induced by the interaction between 
them. 

It is well established that in amorphous solids TLS are coupled with 
phonons. In dielectric glasses at T <C &d (@d is the Debye temperature) 
coupling with acoustic phonons plays a main role. The virtual exchange of 
TLS by acoustic phonons results in the indirect interaction between TLS, 
which decays with the distance R as R~ 3 . In particular, this interaction 
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gives rise to the flip-flop transitions between two TLS (see Fig. 1). The 
transition amplitudes also decay with the distance as it! -3 . 7 ' 8 In a number 
of dielectric glasses, TLS possess their own dipole moment. In this case the 
electrostatic dipole-dipole interaction can become the dominating interaction 
between TLS. 9 This interaction also decreases with the distance as R~ 3 . It 
was experimentally discovered by Arnold and Hunklinger 10 ( see also 9 ) that 
the 1/E? interaction contributes to the spectral diffusion of the TLS energy. 
This spectral diffusion involves the dynamic fluctuation of the TLS energy 
due to its interaction with neighboring TLS, making transitions between 
their levels. 11 

Along with the interaction significance for the acoustic 10 ' 9 and 
optical 12-14 hole burning experiments, it also affects the decay of coher- 
ent echoes. 15 Also the interaction between TLS has been revealed in the 
non-equilibrium dielectric measurements, 5 ' 16 giving rise to the reduction of 
the spectral density near zero energy. This is interpreted as the dipole gap 
formation. 17 In addition, it has been found that the transition from coherent 
to incoherent tunneling takes place if the typical interaction energy exceeds 
tunneling splitting. 18 Recently discovered anomalous low-temperature be- 
havior of amorphous solids in the external magentic field 19-24 can also be 
associated with the TLS interaction. 

The relaxation in a disordered system can be due to delocalized ex- 
citations. The single particle energy derealization is not efficient for the 
ensemble of interacting TLS because the static energy disordering is very 
strong there. Therefore, the Anderson localization 25 of all excitations takes 
place. In other words, the localization in the disordered TLS system is due to 
a large energy level mismatch for a typical pair of TLS (Fig. 1A) compared 
with the flip-flop transition amplitude for this pair. Suppose that a TLS 
excitation energy can experience large time-dependent fluctuations. These 
fluctuations can reduce the energy mismatch value, stimulating level crossing 
and supporting flip-flop transitions (see the resonant pair in Fig. IB). Such 
fluctuations can be induced by either the external alternating (measuring) 
field or they can be due to the spectral diffusion, caused by transitions of 
neighboring TLS. Thus, such fluctuations can promote the derealization in 
the system. In the absence of the external field this derealization can be 
treated as a self-consistent process in which the spectral diffusion induces 
TLS transitions. In its turn, TLS transitions cause the spectral diffusion 
due to TLS interaction. 

This paper is a comprehensive presentation of our results concerning the 
effect of R~ 3 interaction on the dynamic properties of amorphous solids. It 
is organized as follows. In Sect. 2 we address the single particle localization 
problem for zero temperature and zero external alternating field. In Sect. 3 



The Long Range Interaction and the Relaxation in Glasses at Low Temperatures 

the effect of an external alternating field on the single particle localization 
is studied. We demonstrate that the field having a sufficiently large ampli- 
tude and a small frequency can break the Anderson localization. In Sect. 4 
we study the finite temperature many-body delocalization of a TLS energy 
within the TLS ensemble induced by the self- consistent dynamics of res- 
onant pairs (see Fig. IB). Sect. 5 is devoted to the most complicated, but 
experimentally relevant regime where both external field and interaction of 
TLS must be taken into account simultaneously. In summary we present the 
brief discussion of the application of our theory to the experimental studies. 



2. Anderson localization of the TLS energy 

Following Ref., 1 we accept that the distribution function for the bias 
asymmetries A and tunneling amplitudes Ao (see Fig. 1) obeys the universal 
distribution 

P(A,A ) = -£-. (1) 
An isolated TLS can be described by the standard pseudospin Hamiltonian 

hi = -AiSf - A 0i Sf. (2) 
The interaction between TLS can be expressed as 9 ' 30 

V = l'£U{R ij )S!S], U(R ij ) = -^ r , (3) 

where Rij is the distance between two TLS and U is the characteristic 
coupling constant. 

Consider a relaxation of a sole excited TLS with certain energy split- 
ting E = (A 2 + Aq)2. The relaxation channel independent of phonons is 
hopping of an excitation from excited TLS to another TLS with parameters 
A', A' ,E' = (A' 2 + Aq 2 )2, which is initially in its ground state (see Fig. 
IB). As a result, first TLS goes into the ground state, while second TLS 
becomes excited. The inverse process also takes place so that the pair of 
TLS under consideration can be detected in one of the two states separated 
by the small energy | E — E' \ . In what follows, these two states of a TLS 
pair will be referred to as a flip-flop configuration. Such a TLS pair can 
be considered as a new type of the two-level system with the asymmetry 
A p =| E — E' |. In addition, as follows from Eqs. (2), (3), the tunneling 
amplitude coupling two states of the flip-flop configuration of a TLS pair is 
given by the relation 8 ' 26 ' 27 



(4) 
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Thus one can describe the evolution of initially excited TLS in terms 
of the transition between the two states of the flip-flop configuration of a 
TLS pair (Fig. IB). If a certain two- level system with the parameters A 
and Ao was initially, say, in the first state, the time-dependent probability 
W-z{t) to find the system in the second state is given by the solution of 
the corresponding time-dependent Schrodinger equation (hereafter we set 

n = k B = i) 

W2{t) = A^TAl sin V( A2 + A o))- (5) 

This relation applied to the flip-flop transition within the TLS pair (see Fig. 
1) means that the transition efficiency is noticeable provided the parameters 
of a TLS pair obeying the resonance condition 

A p =| E-E' |< Aop(R), (6) 

where the pair transition amplitude Ao p is defined by Eq. (4). Hereafter 
such a TLS pair is referred to as a resonant pair(RP). Then, the frequency 
of quantum mechanical oscillations of population between the potential wells 
of the flip-flop configuration given by (cf. Fig. IB, Eq. (5)) 

r^ 1 = ^(A 2 + A 2 p ) « A 0p (R) (7) 

is completely determined by the pair flip-flop transition amplitude Ao p (i2) 
(4). The inverse of this frequency defines the transition time r. 

The Anderson derealization and/or localization are explicitly associ- 
ated with a TLS capability to form resonant pairs. If the typical TLS has 
on average the large number of resonant neighbors, one would expect the 
derealization. In fact, each resonance shares the energy between two TLS. 
One can use a percolation analogy, connecting all resonant pairs into a clus- 
ter. When the infinite cluster can be formed, excitations become delocalized 
(cf. Ref. 28 ). It is clear that we need to have more than one resonant neigh- 
bour for each TLS (2-3 should determine the threshold 28 ) to make such a 
cluster. When the average number N of resonant pairs formed by one TLS 
is much less than unity, this TLS cannot share its energy with anybody (Fig. 
1A) or it has one resonant neighbour (Fig. IB), with whom it hybridyzes 
its excitation. The probability to add more resonances decreases by a small 
factor of N <C 1, so all excitations are localized. 

How many resonances does a given excited TLS with the parameters 
Aq,E form with other TLS! To answer this question, one should calculate 
the average number N[Aq,E] of TLS, forming a resonant pair with selected 
excited TLS. Using Eqs. (1), (4), (6), we find 



N[A ,E] 
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} 



B. Resonant pair 




Fig. 1. Energy transport between excited TLS (left) and nonexcited TLS 
(right). It is not efficient in the off-resonant case (A), where the energy level 
mismatch exceeds the resonant interval indicated by the figure bracket. For 
the resonant pair (B) the interaction of TLS induces flip-flop transitions 
mixing up two possible quasi-degenerate states, where either first or second 
TLS is excited, while the remaining TLS is in its ground state. 
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r • 0(A Op (r) - A p )0(r - a)0(L - r) 




(8) 



Here © is the Heaviside function involving the resonant condition Eq.(6) 
and minimum and maximum distances between interacting TLS (®(x) = 1, 
x > 1, Q(x) = 0, x < 0), L is the size of the whole system and a is 
the minimum distance between two TLS. Since in all glasses the parameter 
X oc PU is very small (see e. g. 29 > 30 ) 



the number of resonant neighbors is much less than unity for any reasonable 
sample size L. Therefore only a small fraction of two level systems belongs 
to resonant pairs, while majority of the others are immobile. This proves 
the full Anderson localization of all excitations in this system. 



Consider the effect of an external alternating field on the Anderson 
localization within the TLS ensemble. This problem is experimentally im- 
portant because almost all mesuring techniques use the external alternating 
field to probe the system response. At ultralow temperatures T < 0.1K it is 
very difficult to make the field small enough in order to avoid its nonlinear 
effect. It turns out that the interaction-induced relaxation is also affected by 
the field. In this subsection we analyze the simplified single-particle problem 
and a more realistic case will be treated later in Sect. 5 of this paper. 

An approximate zero temperature condition means that one can neglect 
the interaction between excited TLS because their total number is negligibly 
small in this case. In particular, there is no spectral diffusion. Consider the 
effect of the alternating field on the energy spectrum of excitations. 31 Due to 
the field, the energy splitting A; acquires the oscillating part (iiCos{u)t). Here 
a,i is a coupling energy between TLS in the site i and the external field, that 
is a product of the field strength (electric field E or elastic stress e) and the 
coupling constant (TLS dipole moment fi or strain tensor 7, respectively 9 ). 
We assume that the field varies sufficiently slow and the amplitude ctj is small 



compared to the typical scale of the excitation energy E^ = v/A? + Aq, 



X = 4ttPU < 1(T 3 - 1(T 2 , 



(9) 



3. Derealization of Floquet states 




co <C a« <C Ei 



(10) 
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The assumption of a small field amplitude is needed in order to treat the 
field as a weak perturbation. The assumption of low frequency is satis- 
fied in the vast majority of acoustic and dielectric experiments in glasses, 
because the minimum field amplitude usually exceeds 10~ 4 K while the maxi- 
mum frequency is usually below 100kHz corresponding to the energy 10~ 5 K. 
Moreover, usually a ratio a/uo exceeds unity by several orders of magnitude. 

Let bi be the amplitude of the excited TLS state at the site i. An 
evolution of a single-particle excitation can be described by the Schrodinger 
equation with the external field involved 

ibi = (Ei - ai cos (ut)) h + Uijbj. (11) 

3 

Here Uij are hopping amplitudes (4) between pairs of sites i and j. The sub- 
system of "active" or resonant TLS effectively responsible for the dynamics 
consists of TLS with Aj < A(jj. For this subsystem, the factor ^ ^ 1S 
close to unity and can be omitted. Then, considering only the subsystem of 

"active" TLS, we put , , , , 

\U ij \ = U/Rfj (12) 

and assume the uniform density for the energy levels P (E) = P. 30 

Equation (11) can be rewritten in terms of the Floquet state amplitude 
di n (see Appendix) 

zdin = (Ei nu>) din ~\~ ^ ' Ti n -jpdjp, Ti n -jp = UijJp—n ( j (13) 

jp 

with e being an eigenvalue of a Floquet eigenstate quasi-energy. Here J p (x) 
is the Bessel function of order p. 

Equation (13) resembles those usually considered in studying the Ander- 
son localization when the disorder is static. If the hopping term Ti n] j p is ne- 
glected, the eigenstates of the system in the external field correspond to local- 
ized single-site excitations of a TLS coherently dressed by n " quanta" of the 
frequency wits position and the number of the dressing quanta. Therefore, 
the problem is reduced to investigating the derealization of these dressed 
excitations (DE). 

Suppose that an excitation is created at the site i. It can be described 
as a DE with certain double indices iO. To leave the state iO for another DE 
state jn, the condition of resonant coupling should be satisfied, as shown in 
Fig. 2, 



Uij J n 



ai 



(14) 



\Ei — nuj — Ej\ < Tjo-jn 

This criterion is similar to multiphoton resonances in the nonlinear optics. 
The successful transition between two states can be accompanied by an 



Alexander L. Burin and Il'y a Ya. Polishchuk 



A. No External Field 




B. With External Field 



■ n=2 




Fig. 2. The formation of a resonant pair from the initially non-resonant one 
(A) in the external alternating field. The coherent excitation dressing (B) 
enables transitions with the energy change given by the integer number n 
of field frequences. The case n = — 2 satisfies the resonance conditions ( see 
Eq. (14)). 



The Long Range Interaction and the Relaxation in Glasses at Low Temperatures 



absorption or an emission of n quanta, to provide the energy conservation. 
These processes are efficient, when the field amplitude is large to ensure the 
strong nopnlinearity like in our case. 

The local amplitudes a\ vary from site to site due to fluctuations in 
coupling constants. Let us denote the average of a difference |a« — a,j\ by a 
having the same order of magnitude as amplitudes aj. 

The argument of the Bessel function in Eq. (13) is a large parameter 
of order of a/u> » 1 (see Eq. (10)). If \n\ > a/w, the magnitude of the 
corresponding Bessel function is exponentially small, entailing a negligible 
probability of resonant coupling. The opposite condition a >| n \ uj means 
that, in order to have resonant coupling, the energy difference should be at 
least less than the field amplitude a: 

| n |< a/uj — ► \Ei — Ej\ < a. (15) 

On the classical language this condition means that the interaction of TLS 
can be efficient only when the external field stimulates their real level cross- 
ing. 

Under condition (15) and a/uo 3> 1 the Bessel function in (13) can 
be approximately replaced by its asymptotic value. Omitting the standard 
oscillating prefactor which plays no role in a random discrete problem, one 
can estimate the coupling amplitude for DE Eq. (12) as 

(16) 

Consider the derealization of an excitation due to the alternating field 
within the framework of the concept of resonant coupling. Two sites i and 
j are in the resonance when the condition Eq. (14) is satisfied at least 
for some n < a/ui (see Eq. (15) and Fig. 2). Since the energy splitting 
between two subsequent levels is equal to the field frequency u, the resonance 
always occurs if the hopping amplitude T^j n exceeds uj and condition (15) 
is fulfilled. As directly follows from Eq. (16 ), this happens when sites i and 
j are separated by a sufficiently small distance Rij 

Rij < n = {U/V^) 1/3 - (17) 

Thus all TLS j located closer than the crossover distance r* to TLS i under 
consideration with \Ei — Ej\ < a are resonantly coupled with each other. 
Taking into account that the number of such centers in a unit volume is 
equal to Pa, we find that the total number of the resonant neighbors within 
the sphere of a radius r* is given by 



W (r*) = 
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At larger distances R > r* the coupling amplitude (16) is less than to. Then 
one should treat separately all a/uj possible resonances (see Fig. 2). Thus 
for the given pair the probability of resonance increases by the factor of a/uj 
due to the increase in the number of possible resonances. On the other hand, 
the probability of each resonance decreases by the factor y/uJJa because of 
the reduction in the flip- flop transition amplitude (see Eq. (16)) compared 
to the zero-field case Eq. (4). As a result, the total number of resonances in 
each layer r* < r\ < R < 

X* = 4itPU x [^, (19) 
V u> 

increases by the factor y 1 ajuo as compared to the zero field case (8). Thus, 
in the precence of the external field the parameter W (r±,r2) increases by 
the large factor yja/uj. 

According to the Levitov's analysis of the derealization problem in the 
case of the 1/ R? energy hopping amplitude, 32 the parameter x* is a decisive 
parameter for the derealization. As follows from Eqs. (18) and (19), when 
X* <C 1 the first resonance occurs at the distance R a ~ r^e 1 /** 3> r*. Thus, 
the first hop of the excitation occurs to the distance R a , taking some time 
t\ determined by the inverse hopping amplitude t\ ~ R a /U oc e 3 /**. Then, 
the time required for the second hop, when the next resonance appears, 
is exponentially large compared with t\? 2 Therefore, the derealization is 
exponentially slow, if any. 

Consider the opposite case %* > 1. Then, the number of resonances 
exceeds unity in each spherical layer n < R < 2r± (n > r*) and the reso- 
nant sites form an infinite cluster meaning an existence of the delocalized 
state. 8 Thus, when a ratio of the external field amplitude to the frequency 
is sufficiently large to provide the condition 

X* > 1, (20) 

the derealization of excitations takes place. 

The inverse time of a single hop between two nearest resonant neighbors 
in the derealization regime can be treated as a relaxation rate for the DE 
located at a certain site. The distance R a < r* between these neighbors can 
be estimated from the relation W(R a ) w 1 (see Eq. (18)). This radius is 
determined as 

R a « (Pa)- 1 / 3 (21) 

and the typical dipole-dipole hopping amplitude corresponding to the dis- 
tance R a is given by Eq. (16) 
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T(R a ) = ^^/~anu X *. (22) 

"■a 

This expression gives an estimate of the inverse lifetime (the relaxation rate) 
for the BE at an arbitrary site 

t* 1 ~ UJX* = {auj) 1/2 X (23) 

4. Many body derealization, dephasing and relaxation 

Here we will discuss the case of the zero field and a finite temperature 
T > 0. The Anderson localization of excitations proved in Sect. 2 takes 
place only if every resonant pair RP can be treated independently. This 
assumption is valid provided that one can neglect the interaction between 
different excited TLS. We will see that this is not the case at any finite 
temperature T > 0. Below we will show that at any finite temperature the 
long-range 1/R 3 interaction of excited TLS leads to the irreversible dynamics 
and relaxation. This relaxation is essentially of a many-body origin. In other 
words, one should take into account simultaneous transitions in two or more 
RP (see Fig. 3) and thus at least four TLS will participate in an elementary 
process. 8 ' 27 ' 33,34 

Any RP has four energy levels. Two of them correspond to the flip- 
flop configuration mentioned above (see Fig. 2B). The two other states of 
a pair correspond to the configuration where both TLS are either in their 
excited or their ground states. In fact, the typical energy of excited TLS 
in a resonant pair is given by the thermal energy T, while their flip-flop 
transition amplitude Eq. (4) is generally much smaller due to the weakness 
of the interaction Eq. (8). Therefore the flip-flop interaction Ao p can connect 
only pairs of TLS, where one of them is in its excited state and the other 
one is in its ground state (Fig. 2). The other two states are separated from 
the flip-flop pair by the large energy gap of order of the temperature. 

The states of the flip-flop configuration are separated by the energy 
interval A p =| E—E' |. In spite of the fact that E, E' ~ T, one can construct 
RP with A p ,Aop <C T. Then, even if the interaction V{R) between these 
TLS is weak, the condition A p < Aq p can be valid. Therefore the two levels 
in the resonant pair can be strongly coupled (see Eq.(6)). In the remainder 
of this paper we will consider only flip-flop configurations of RP. We will 
treat RP as a new kind of the two-level system with the energy asymmetry 
A p =| E — E' | and the tunneling amplitude Aq p {R) (4). 

Thus below we investigate the relaxation of this novel RP type of the 
two-level system for which the distribution function of parameters A p and 
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Ao p is defined as 8 ' 30 

P( 2 )(A P , A 0p ) = (S(A p - | E-E> |)«5(A 0p - |>^)>- (24) 

The brackets denote two averaging, namely, thermal averaging and averag- 
ing over the distribution of parameters of original TLS (see Eq. (1)). In 
addition, the integration over the distance R includes all possible pairs, so 
that the left-hand side of Eq.(24) can be rewritten as 

J d 3 R J P(A, A )dAdA J P{A',A' )dA'dA' n(E)- (25) 

(1 - n(E'))5(A p - | E-Ef \)S (a 0p - , (26) 

with n(E) = [1 + expiE/T)}- 1 being the probability to find TLS in its 
excited state. To evaluate Eq. (26), we take into account that the integral 
is determined mainly by E ~ E ~ T. Therefore, due to A p <C E ~ T 
one can omit A p in the argument of the first delta-function. Then one can 
easily evaluate the pair distribution function Eq. (26) within the logarithmic 
accuracy 8 ' 33 

P^(A p ,A 0p )^(PT)(PU)^. (27) 

Note that, the (Ao p )~ 2 dependence in the last expression is a consequence 
of the i? -3 interaction and it is not related to the initial distribution Eq. 
(1). Although the distribution function Eq. (27) differs from Eq. (1), it 
also does not depend on an asymmetry parameter A p . The pair distribution 
function Eq. (27) has a stronger singularity at small Ao p than the original 
TLS distribution function Eq.(l) at small Ao . For this reason, the con- 
centration of low energy RP excitations is larger (and, correspondingly, the 
average distance between them is smaller) than for initial TLS. The coupling 
constants U(R) both for TLS and RP are of the same order of magnitude, 
so that low energy RP interact stronger than TLS (See Fig. 3). The reason 
for that is that having iV 3> 1 TLS we can make N 2 3> N pairs of them. 

The ensemble of RP, new kind of TLS, is described by the initial Hamil- 
tonian Eqs. (2), (3). The only but a key distinction from the initial model is 
that the distribution function (27) should be used instead of Eq. (1) at finite 
temperature T > 0. Remember that resonant pairs are only those pairs for 
which A p < Aq p . Therefore each resonant pair has the single characteristic 
energy given by its flip-flop amplitude Ao p . Since resonant pairs are mainly 
formed by TLS with A ~ Ao ~ T, the charactreristic pair transition ampli- 
tude can be estimated as U/R 3 (4), where R ~ (U/Aop) 1 / 3 is the size of the 
pair. 
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A. Initial State 




C. Forbidden Final State 




Fig. 3. The interaction of resonant pairs leads to the many-body dereal- 
ization (A). The excitation is transferred due to the simultaneous flip-flop 
transitions of four TLS to their new states (B). The transitions involving 
two TLS belonging to different resonant pairs are not allowed because of the 
energy mismatch (C). 



Resonant pairs interact with each other. We will show following 
Refs. 8 ' 33 that this interaction inevitably leads to the energy derealization 
induced by collective flip-flop transitions similar to that shown in Fig. 3. 
It is convenient to separate the whole sequence of resonant pairs into the 
infinite set of strips k = 1,2,3, ... These strips are defined in the space of 
characteristic energies of resonant pair. Each strip k is formed by resonant 
pairs having transition amplitudes within the range (Ao p (k) /2, Ao p (k)), with 
Aop(^) = T/2 k . One can estimate the concentration of RP within the 
strip k making use of the distribution function Eq. (27) 

N k « (PT){PU) = N*. (28) 

It is important that the RP concentration in each strip k is completely 
independent of the value Ao p (fc). Separating strips corresponding to all 
permissible values Aq p , one can cover completely the whole ensemble of 
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resonant pairs. Within any strip the RP concentration is constant iV* Eq. 

(28). Therefore the average distance between RP within any strip does not 

_ i 

depend on the kind of it. This distance is equal to R* ~ iV* 3 . Since the 
interaction between RP is of the same origin as the interaction between TLS, 
it is given by the same expression U/Rf (see Eq.(3)). Therefore, the typical 
energy of the interaction between RP within any strip is given by 

£/* « U/Rl « T X 2 . (29) 

Since RP is a kind of TLS, one can introduce a concept of a flip-flop 
configuration for two RP shown in Fig. 3 and derive an expression for the 
transition amplitude Aq 4 "* between the levels of this configuration (which is 
similar to that of Eq.(4)). Because interacting pairs are both resonant this 
transition amplitude is given by their characteristic interaction Eq. (29) (cf. 
Eq. (4)). Consider the effect of this flip-flop interaction on RP belonging to 
different strips. 

We should start with fastest, high energy resonant pairs of the strip 
k = 0, while other pairs are slower and can be treated quasi-statically. It 
is clear that their characteristic interaction Eq. (29) is much smaller than 
their characteristic energy T because of the weakness of the TLS interaction 
X < 1 Eq. (9). Thus, they are isolated from each other. Their dynamics 
is faster than for other pairs having smaller energies, so they can be treated 
adiabatically with respect to other RP. Then, the fast RP possibly contribute 
to a small redefinition of parameters for slower pairs without any effect on 
their average properties. Then to analyze the interaction effect for other pairs 
we should go all the way down in energy, increasing a strip index k. It is 
clear that the interaction becomes significant when the characteristic energy 
of pairs Aq/c approaches the interaction energy Eq.(29). Remember that the 
interaction of RP belonging to each strip is given by the universal value 
Eq. (29). The pairs with energy Aqa; < U* cannot be treated independently 
because their interaction is greater or equal to their characteristic energy. 

The strips with the energy smaller than [/* satisfy the conditions of 
the Anderson derealization because the characteristic flip-flop transition 
amplitude £/* between RP is greater or equal to their energy disordering, 
determined by their characteristic energy Aofc- The elementary process of 
an energy hop between pairs is demonstrated in Fig. 3. It is important that 
in order to transfer the energy, the simultaneous transition of all four TLS 
is required, while the transition of two TLS belonging to different resonant 
pairs (Fig. 3C) is forbidden because of the energy mismatch of order of the 
thermal energy T. Thus RP excitations in the cluster formed by the strip 
with the energy ?7* are delocalized. The relaxation rate of excitations within 
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this strip is given by the inverse characteristic interaction 

r" 1 ~ U* « T{PU) 2 . (30) 

It turns out that the time r* defines the dephasing time within the 
ensemble of TLS. Consider the time-dependent fluctuation <53>(i) of the phase 
at arbitrary probed TLS, induced by the energy-splitting fluctuation 5E(t). 
The proper dephasing time t 2 can be determined from the relationship 

5<5>{t 2 ) = t 2 5E(t 2 ) « 1. (31) 

To show that the dephasing time is given by the time r* , let us calculate the 
energy-splitting fluctuation 6E(t) within the time interval t = r*. The overall 
density of RP and, therefore, approximately the density of TLS undergoing 
the transition during the time t = t* is of order of the RP density N*. 
Bearing in mind that every RP experiencing the transition at the distance 
r from probed TLS contributes a value U(r) (see Eq. (3)) into the energy- 
splitting fluctuation for probed TLS, one can estimate 

8E(n) = U(r) « N*U « 17(12.). (32) 

r 

This estimate was made taking into account that the sum in Eq. (32) is 
defined by the largest term coming from the closest neighbour contribution. 11 
THis nearest neighbour is separated from TLS under consideration by the 

—1/3 

distance R* = N* 

Then, taking into account Eq. (30), one can find that t~ 1 5E(t*) 1. 
The comparison of this relation with Eq. (31) allows us to link r* with t 2 
and therefore 

r 2 l « U{R*) « T X 2 (33) 

This dephasing rate decreases linearly with the temperature. Compared with 
the phonon- induced channel providing a T 2 -dependence (see Ref. 11 ), the 
dephasing governed by the mechanism concerned predominates at sufficiently 
low temperatures. 

Regardless of the relaxation mechanism, the spectral diffusion rate is 
defined as 9,11 



E 



_ -2 



(34) 



On the other hand, in a system with the 1/i? 3 interaction this rate can be 
expressed as the ratio of the spectral diffusion amplitude %T and thermal 
TLS relaxation time Ti 9,11 



rf = X Tlr x . 



(35) 
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Eq. (35) is based on the time evolution of a TLS energy perturbed by 
transitions of neighboring TLS in the system with \/R? interaction. If the 
characteristic relaxation time of thermal TLS having the density PT is given 
by ri, then during the time t < t\ the fraction of all TLS rir(i) ~ PTt/ri 
undergoes transitions. The characteristic perturbation in each TLS energy 
caused by these transitions is defined as Unx{t) = UPTt/r\. The phase 
fluctuation can be estimated using Eq. (31) 5$ ~ UPTt 2 /n. Assuming 
that at t = T2 the phase fluctuation becomes of order of unity (see Eq. (31)) 
we derive Eq. (35). This derivation is correct, when TLS transitions have 
no correlations with each other. This is of course true for phonon stimulated 
relaxation, while in our case the spectral diffusion should be treated self- 
consistently. It was shown in the previous works 8,33 that despite of certain 
memory effects in the spectral diffusion Eq. (35) remains approximately 
valid. 

Using (33) and (35) one can estimate the relaxation rate for thermal 
TLS as 8 

rf 1 = X 3 T. (36) 

5. Relaxation at finite temperature and strong field 

Previously, we have shown that the 1/-R 3 interaction of excitations and 
the external alternating field can independently result in the derealization 
and in the irreversible TLS dynamics in the regime of arbitrary strong dis- 
ordering. Under conditions of the real experiment both the TLS interaction 
and the external alternating (measuring) field are significant. Below we con- 
sider the combined contribution of two those effects on the TLS relaxation. 

Our analysis starts with the very low frequency of the external field, 
where the field can be treated as static. Then the frequency will substantially 
increase to its maximum value u> ~ T. The case of larger frequences u> > T 
is beyond the scope of this paper because of almost unavoidable heating. 
The results of our consideration are placed into the Table 1. 



Table 1. Frequency-dependent relaxation rate. 



Frequency range 


rf 1 




0<u;< {T X 2 Y /a 


T X 3 


Tx 2 


(T X 2 ) 2 /a < oo < a X 2 


(aio) 1 x 


(«^T 2 ) 1/4 X 


ax 2 < u < Tx 


ax 2 


(aT) l / 2 X 3/2 


Tx<to <a 


(a/io) T X 3 


{a/u)V 2 T X 2 


a <co <T or a < Tx 




Tx 2 



Below we consider the case of a relatively large field amplitude a > Tx 
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where the field effect is significant. One can show, using similar methods 
as below, that for smaller field amplitudes a < T\ the field effect can be 
ignored and the relaxation is defined by the results of the previous section 
Eqs. (33), (36) for the relaxation and dephasing rates, respectively. It is 
interesting that the similar non-linearity criterion has been recently found 
in metallic glasses. 35 

Remember, that the parameter T\ = U (PT) determines the character- 
istic interaction of thermal TLS (Aq ~ A ~ T). It represents the amplitude 
of the interaction induced spectral diffusion. The condition a > T\ means 
that the fluctuation of a TLS energy induced by the external field exceeds 
the interaction-induced fluctuation. We also assume that the field amplitude 
is less than the thermal energy to avoid heating. Thus in our consideration 
below the field amplitude a satisfies the inequality 

T X <a<T. (37) 

5.1. Quasistatic field 

When the frequency of the external field u is very small, the field is 
almost static and we can ignore it. Then the relaxation and dephasing 
rates are defined by the zero field limit Eqs. (33), (36), respectively. In 
this regime the fluctuations of a TLS energy are due to the self-consistent 
spectral diffusion alone. The spectral diffusion rate is given by Eq. (34). 
The rate of a TLS energy fluctuation induced by the external field is coo. 
The external field can be treated as static until the spectral diffusion rate 
exceeds the fluctuation rate due to the field tvT 2 > auj, which takes place at 
^ < (Tx 2 ) /a (see the first row in the Table 1). 

5.2. Adiabatic relaxation induced by slow field 

At larger frequency 

u > (T X 2 ) 2 /a (38) 

the field cannot be ignored. It can stimulate transitions of TLS and irre- 
versible relaxation due to energy level crossings of different TLS. At a large 
field amplitude Eq.(37) crossing between two TLS energy levels Eq and Ej 
( Eq — Ej — a cos ut = 0) 36 is possible provided that 

\E -Ej\<a. (39) 

The average distance between two TLS satisfying Eq. (39) is given by R a ~ 
{Pa)~ 1/3 (see Eq. (21)). 

Consider a flip-flop pair of TLS with the close energies satisfying Eq. 
(39) and separated by the distance Rj. This pair undergoes level crossing 



Alexander L. Burin and Il'y a Ya. Polishchuk 



during the field oscillation period u 1 . If the pair transition amplitude is 
large V 2 = (u/Rfj > E ~ au>, the adiabatic relaxation takes place. Then 
according to the Landau - Zener theory the excitation will transfer within 
a TLS pair with the almost unity probability if the size Rj of a pair is 
sufficiently small (see Eq. (17)) 

Rj < tv ~ (U/^) 1/3 . (40) 

The adiabatic regime takes place when each thermal TLS has the large 
number of adiabatic neighbors, satisfying Eqs. (39), (40). This requires 
?V > R a , which leads to the frequency constraint from top (see Table 1, the 
second row) 

to < a X 2 (41) 

For given thermal TLS the relaxation rate r-f 1 can be estimated as a fre- 
quency of adiabatic level crossings given by the product of the field oscillation 
frequency to and the number of adiabatic level crossing events during the sin- 
gle field oscillation. The latter number is given by the number of adiabatic 
neighbors x*, satisfying Eq. (39) and located within the sphere of the radius 
?V around given TLS, x* — Par^ = (a/w) 1 ^ 2 x » 1 

r 1 - 1 ~ X ^ = (aa;) 1/2 x (42) 

It is interesting that the relaxation rate (42) coincides with that of Eq. (23). 

In the adiabatic regime we can estimate the dephasing rate using Eqs. 
(35), (42) 

Tz 1 = (Tx/n) 1 ' 2 = (tuJt 2 ) 1 ' 4 X (43) 

It is important that dephasing (43) is so fast that the phase coherency be- 
tween the periodic events of level crossings occurring during the period uj^ 1 
can be ignored. This is true in the whole adiabatic frequency domain (the 
second row of the Table 1) because the dephasing time is in fact less than 
the oscillation period (see Eqs. (37), (41), (43)) 

ojt 2 = (^/(ax 2 )) 3/4 (a/T) 1/4 x 1/2 < 1. (44) 

5.3. Non-adiabatic regime at intermediate frequencies 

Consider the nonadiabatic regime, which takes place at higher frequen- 
cies 

lo > a X 2 . (45) 

This condition is opposite to Eq. (41). In this regime the vast majority of 
thermal TLS undergo nonadiabatic level crossings (see Eq. (46) below). 
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In the adiabatic case we have dealt with the regime 00T2 < 1 (see Eq. 
(44)). In the nonadiabatic regime characterized by Eq.(45) one should dis- 
tinguish between two cases LOT2 < 1 and UT2 > 1 where the phase memory 
between two subsequent level crossings either exists or does not exist, re- 
spectively. 

Consider the first case LOT2 < 1 valid at the border of adiabatic and 
non- adiabatic domains u> = ax 2 (see Table 1). Here dephasing is fast and 
phase correlations between periodic level crossings can be ignored. Most 
efficient flip-flop transitions occur between thermal TLS separated by the 
average distance R a Eq. (21). In fact this distance is the characteristic 
separation of neighboring two-level systems with energy difference less than 
the field amplitude a, required for level crossing. On the other hand at larger 
distances a nonadiabatic transition probability decreases with the distance 
very fast as V 2 ~ i?~ 6 . Therefore the most efficient energy transfer should 
occur between nearest neighbors. The relevant transition amplitude is given 
by Eq. (12) V a = U/R^ ~ a\- In the large frequency case Eq. (45) one has 

V 2 < E = auj (46) 



and transitions are nonadiabatic. According to the Landau - Zener theory, 
their probability per one energy level crossing, occurring a few times for the 
period w -1 , is W a ~ V2I (aw) = (a/a;) x 2 ■ This probability defines the TLS 
relaxation rate as the inverse average time between two successful transitions 

rf 1 = W a u ~ a X 2 . (47) 



Then, using Eq. (35), one finds 

r, 1 = a^T^x^ 2 (48) 



The above derivation remains valid until breaking the condition Eq. (44) at 

u > a i/2ri/y/2. (49) 



In the opposite case of L0T2 > 1, the field periodicity is significant so 
that the formalism of dressed excitations (DE) introduced in Sect. 3 (see 
also Ref. 31 ) becomes applicable. 

The further analysis depends on the relationship between DE energy 
splitting uj and the spectral diffusion amplitude Tx (see Eq. (37)). We start 
with the case of lower frequencies (see Table 1) 



uj<T X - 



(50) 
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where all pairs of TLS with an energy difference less than the field amplitude 
a inevitably undergo DE level crossing during the spectral diffusion quasi- 
period r\. These crossing levels are coupled by the transition amplitude of 
dressed excitations Eq.(16). 

In the regime of Eqs. (49), (50) the relaxation is induced by non- 
adiabatic level crossing caused by the spectral diffusion. On the other hand, 
the spectral diffusion is caused by the relaxation dynamics of TLS, so the 
process is self-consistent. 8 ' 27 ' 33 ' 34 

Assume that there is an existing characteristic relaxation rate of thermal 
TLS T-f 1 . Transitions of thermal TLS change other TLS energies giving 
rise to the spectral diffusion. The spectral diffusion leads to level crossing 
of DE, stimulating irreversible transitions with the output rate r out . The 
self-consistent relaxation mechanism requires the input rate r-f 1 to be equal 
to this output rate. 

The rate of transitions induced by the spectral diffusion can be esti- 
mated as the number of transitions induced by level crossings during the 
quasi-period of the spectral diffusion n multiplied by the frequency of spec- 
tral diffusion cycles r-f 1 . Energy level crossing for two DE happens when 
(cf. Eq. (14), Fig. 2) 

E l -E 2 = nu (51) 

(where n is an integer number). The total number of such crossings due to 
the spectral diffusion induced energy fluctuation T\ for the time t\ is given 
by 

N T1 ~ (T X /u). (52) 

Multiple crossings for the same TLS pair give a logarithmic correction to Eq. 
(52) and can be neglected in the qualitative scaling approximation. 8 ' 33,34 
This is due to the specific behavior of the spectral diffusion induced by 
the 1/i? 3 interaction. For this spectral diffusion the characteristic energy 
fluctuation is directly proportional to the time. 9 ' 11 Therefore this is the 
anomalous diffusion process compared to the normal diffusive t 1 / 2 behavior 
of the displacement. For this super- diffusion case the probability of energy 
to return back to its initial value is not so large as in the case of the normal 
diffusion. 

We consider nonadiabatic flip-flop transitions between only neighboring 
thermal DE separated by the distance R a (see Eq. (21)), because the nona- 
diabatic transition probability to larger distances drops with the distance 
very fast (i? -6 ). The probability of the nonadiabatic transition during one 
level crossing for such a pair of thermal DE induced by the spectral diffusion 
with the rate <| i |>~ r^ 2 can be estimated using the Landau-Zener theory 
in the nonadiabatic limit. The transition probability can be found using the 
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transition amplitude Eq. (16) for R = R a 

W* = (V a T 2 ) 2 < 1 

The relaxation rate is expressed as 

rf 1 = r out = N T1 W*/n ~ a X 2 . (54) 

One can show that if we take the input relaxation rate slower than Eq. (54), 
the output rate will be faster than the input one, and, if we take the input 
rate larger than Eq. (54) the output rate will be smaller than the input 
rate. Therefore the solution (54) is stable and the only possible. Thus, the 
nonadiabatic single particle relaxation described by Eqs. (47), (48) takes 
place within the whole frequency domain ax 2 < w < Tx (Table 1, the third 
row). 

5.4. High frequencies 

When the external field frequency u exceeds the scale Tx, the spectral 
diffusion does not necessarily lead to DE energy level crossing for closely 
located pairs, considered in Sect. 5.3. The single-particle derealization in 
this case does not occur similarly to the case studied in Sect. 2 because the 
field frequency is large compared to that needed for the derealization (see 
Table 1, line 2). On the other hand, there exists a finite concentration of 
excited DE (excited TLS dressed by some certain number of the external 
field quanta, see Fig. 2) and many-body relaxation should be similar to that 
described in Sect. 4. 

Following the approach of Sect. 4 one can introduce the concept of the 
resonant pair of DE, instead of RP (resonant pairs of TLS ). A resonant 
pair of DE {RPDE) is defined as a pair of DE separated by the distance 
i?i2 with sublevel energies E\ and E 2 , which obeys the resonant condition 
for some integer n < a/u> (see Fig. 2) 

A pn = |£i - E 2 + nu\ <T(R). (55) 

Here the transition amplitude T (R) ls taken from Eq. (16). Like a 

resonant pair, one can treat RPDE as a new kind of two-level system. The 
parameter A pn is the asymmetry energy for RPDE. Then, strictly following 
the derivation between Eq. (24) and Eq. (28) one can estimate the concen- 
tration of RPDE for any given value of transition amplitude T(R). This 
concentration is given by the expression 



(53) 



N*^T X 2 ^ 

V UJ 



(56) 
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The concentration N* is independent of the given parameter T(R) just like 
the concentration of RP iV* defined by Eq. (28) does not depend on the 
parameter Ao p . The appearance of the factor in Eq. (56) in comparison 
with Eq. (28) can be understood as follows. The concentration N* com- 
pared with N* acquires the factor n* = a/uj since resonance condition (55) 
should be valid for at least one integer \n\ < so that the number of pos- 
sible distinguishable resonances for the single pair multiplies the resonance 
probability by the factor n*. On the other hand, the resonance probability 
decreases by the factor 1 / ^Jnl proportionally to the reduction in the transi- 
tion amplitude Eq. (16). The overall effect is just the increase in the density 
of resonant pairs by the factor \Jajoj (cf. Eq.(19)). Accordingly, the inter- 
action between these resonant DE pairs enhances by the same factor and 
the new dephasing rate is given by this universal interaction 

l/r 2 = T X 2 (a/u) 1 / 2 . (57) 

Making use of Eq.(35) one can estimate the relaxation rate as 

Ti l = (Txriy 1 = (a/u;)Tx 3 . (58) 

When the frequency exceeds the amplitude a < to < T, excitation 
dressing becomes negligibly small and we return to the linear regime (the last 
row of the Table 1). The regime when the frequency exceeds the temperature 
is realized in several echo measurements. 2 This regime can lead to heating 
and requires a special study. 

6. Discussion 

We have described the self-relaxation rate of TLS at different temper- 
atures, external field amplitudes, and frequencies. We found that the relax- 
ation rate is either temperature independent (see the Table 1), or decreases 
with the temperature as T. The phonon stimulated relaxation r-f 1 oc T 3 (see 
Ref. 9 ) is much slower in the low temperature limit and the self-relaxation 
described above should dominate when T — ► 0. 

The suggested theory predicts the parametric dependences of relaxation 
rates in various regimes. In contrast to the theory for the phonon-stimulated 
relaxation (see the review of Hunklinger and Raychaudchary 9 and references 
therein) we are not able to determine the numerical coefficients for each rate 
process shown in the Table 1. This missed factor can be of order of unity, 
but it can also be much greater or much smaller than unity, like 0.01 or 10 2 . 
The example of the latter situation is given by the tunneling rate involving 
two-phonon processes, examined by Kagan and Maksimov 37 for the quantum 
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diffusion problem, where they found the large numerical prefactor of 1000 
due to the large factorial factor involved. 

We do also expect large numerical prefactors for our expressions. There 
are several reasons to have them. First, each factor x possibly involves the 
spherical integration factor Air or in addition to PU factor. In this 

paper we have included the factor 4-7T into the definition of our factor \ Eq. 
(9) contrary to the previous work. This is done because this factor nat- 
urally appears in the localization criterion (8) and correspondingly it can 
enter the rate expressions of the table 1. Thus the parameter of interaction 
weakness PU ~ 10 -4 - 10 -3 can be increased by the order of magnitude. 
Accordingly, the "linear" expression for the relaxation rate Eq. (36) agrees 
better with the results of the systematic internal friction measurements per- 
formed by Classen, Burkert, Enss, and Hunklinger. 6 We believe that the 
factor x = 4nPU reflects the absolute values of relaxation rates better than 
the smaller factor x = PU itself. The accurate calculations of the numer- 
ical factors are beyond the scope of our qualitative study. The alternative 
explanation for a quantitative disagreement 6 can be based on the strong 
non-linearity. For instance, the extra- factor a /to in the relaxation rate (see 
Eq. (58)) can account for the difference of theory and experiment as well. 
The comparison of our predictions with experiment can still be performed 
using the experimental data for different glasses, with different values of the 
parameter x- Note that the preliminary numerical analysis of the Floquet 
state delocalization, described in Sect. 3, in the equivalent one-dimensional 
model supports our expectations of the large numerical prefactor, exceeding 
unity by at least one order of magnitude. 

Another possible problem of the direct application of our theory to the 
experiment is that the interaction-induced relaxation leads to the equilibra- 
tion within the TLS subsystem rather than for the whole system. This is 
similar to spin- spin and spin-lattice relaxations in the NMR problem. Then 
in order to to describe the heat balance between TLS and phonons, one 
should introduce the separate TLS and phonon temperatures and perform 
the thermal balance analysis for the whole system. This study is beyond the 
scope of this paper. 

The crossover temperature between two regimes depends on the exter- 
nal field parameters a and u. It usually reduces with decreasing the field 
amplitude and can also decrease with decreasing the field frequency (see 
the Table 1). This knowledge helps us to understand the absence of the 
interaction-stimulated relaxation reported by Pohl and coworkers in Ref. 38 
In this work, the strain field amplitude was made extremely small e < 10~ 8 , 
while the frequency to ~ 0.5MHz is higher than in other group studies 4-6 
and belongs to the range described by the fourth column of the table 1, when 
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the relaxation rate decreases with increasing u. 

Our results agree with the recent low-temperature measurements of the 
dielectric constant by Ladieu and coworkers, 24,39 which can be interpreted 
assuming that the TLS relaxation rate becomes temperature independent. 
This shows the dependence on the external field amplitude y/a (the sec- 
ond row of the Table 1). Another very interesting discovery of Ladieu and 
coworkers includes the sensitivity of the TLS relaxation rate to the sam- 
ple thickness. The suggested theory is essentially three dimensional. If the 
thickness of the sample will be less than the typical distance between reso- 
nant pairs (around few tens of nanometers) then the irreversible relaxation 
will remarkably slow down. This agrees with the observations of Ref. 24 ' 39 
More accurate theoretical analysis of the data is necessary for the careful 
interpretation. 

Our predictions can be directly verified using the non-equilibrium mea- 
surement technique developed by Osheroff and coworkers, 5 ' 40 which is based 
on the analysis of the system response to the large sweep of the external elec- 
tric field Enc taking different times t s . Our results for the field-stimulated 
relaxation are valid in this case if one takes a ~ Edc^tls (where ^tls ~ 1-D 
is the typical TLS dipole moment) and to ~ 1/t s . Since the maximum pos- 
sible amplitude a is very high a > 0.1 A", while the minimum "frequency" is 
as low as uj ~ Is -1 , all the regimes described in the Table 1 can be attained 
and analyzed. 
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8. Appendix 

To solve Eq. (11), let us introduce the partial amplitudes of the Floquet 

states Ci n so that 

bi = exp [—ist] c in exp (—inuot) (59) 

n 

Substituting this expression into Eq. (11) and taking its Fourier transform, 
we arrive at the equation 

ec im = (Ei - muS) c im - y {c im+ i + c im _i) 

+E c w- ( 6 °) 

3 

Then multiply the last equation by the Bessel function J n - m and per- 
form the summation over m. The result can be presented in the form 

^ Z/m CimJn—m (7^") — (.^i TIU}) X^m ^imJn—m (~) 

+ T, jm UijCjmJn-m (£) + £ m Qm ((n - m) W (61) 

■^n-m (fj) - % (<^(n-m)+l (SO + ^(n-m)+l (SO)) 

Due to the recurrsion properties of the Bessel functions, each term in the last 
sum of the above equation vanishes. Finally, using the Graph summation 
formulae 

Jn—m ( ] / Jp—n ( ) Jp—m ( ] • (62) 

\lo J \uj ujJ \ujJ 

in the sum J2jm i n Eq- (61) and introducing the notation 

din = ^ ' CjmJn—m 
m 

we obtain 

zdin = (ej nui) di n + ^ ' Ui n -jpdjp, 

jp 

U- ■ -U--J ( aj ~ ai ^ (6A) 

V LU / 
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